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We study the composite fermion construction at and below the single vortex {L = A'') state of 
weakly interacting rotating Bose gases, presenting a new method for handling the large number of 
derivatives typically occurring via the Slater determinant. Remarkably, the CF wave function at 
L = N becomes asymptotically exact in the thermodynamic limit, even though this construction is 
not, a priori, expected to work in the low angular momentum regime. This implies an interesting 
mathematical identity which may be useful in other contexts. 
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INTRODUCTION 

Rotating atomic Bose condensates have been subject 
to ever increasing theoretical and experimental interest 
in recent yearsfj. On the experimental side, substantial 
progress has been made in increasing the rotation such as 
to produce ever larger amounts of vortices and, 
very recently, coming close to the point where the vortex 
lattice is expected to melt and the system is expected 
to enter the quantum Hall regime 0. Theoretical stud- 
ies have addressed the entire range of angular momenta, 
from few-vortex states {L ^ N where L is the total an- 
gular momentum and N the number of particles), up to 
the regime L ~ iV^ where the system displays a number 
of close analogies to the fractional quantum Hall effect 
(FQHE)0,I1I3. 

A theoretical tool which has been very successful in 
describing the physics in the quantum Hall regime is the 
composite fermion (CF) construction lOj. Originally de- 
veloped in the context of the fractional quantum Hall 
effect, CF wave functions have, more recently, been ap- 
plied to the high angular momentum states of weakly 
interacting Bose gases0, 0, 0, By construction, 

this method has been assumed to be applicable only in 
the FQH regime, L ~ N'^. However, as was pointed out 
in there were some indications that the CF wave 
functions might do surprisingly well even at the lowest 
angular momenta. In this work we thus discuss the CF 
construction in the regime < L < iV. In particular, 
we present a numerical study which shows that the CF 
wave function for the single vortex (L = N) becomes ex- 
act in the large limit. One of our main motivations 
for studying this case in more detail were earlier calcula- 
tions by Cooper et al [oj, which revealed an increase in 
the overlap as function of N for up to 10 particles. We 
confirm that this trend continues in a systematic way for 
higher A'^, with the deviation of the overlap from unity 
decreasing as ~ Since the exact ground state wave 

function is known analytically in this case 0, lisj . 
this result suggests a non-trivial mathematical identity, 
valid asymptotically in the limit — > cxd, which may 
turn out useful in other contexts. Our solution involves 



a new way of handling the large number of derivatives 
typically occurring in the CF construction when going 
to low angular momenta. This large number of deriva- 
tives has been one of the main obstacles to taking the 
CF construction to higher particle numbers than, say, 10 
[l6j |. We have computed the overlap of the single vor- 
tex CF wave function with the exact one for up to 43 
particles, which is a far higher number than could be 
achieved previously. This was done by taking advantage 
of the mathematical properties of the determinants in- 
volved, which enabled us to rewrite the wave functions 
in a simpler way, without derivatives, before performing 
the numerical calculation of the overlaps. A straightfor- 
ward modification of this method allowed us to study 
other angular momenta as well, and revealed some pe- 
culiar mathematical properties of the CF construction. 
We would like to emphasize that this is mainly a techni- 
cal work in the sense that our main emphasis is on this 
new method and on the mathematical properties of the 
CF construction. While our study was performed in an 
angular momentum regime where exact analytical wave 
functions are known, we hope that our results and meth- 
ods may turn out useful to study states beyond the single 
vortex, where this is not the case. For example, the same 
approach may be used to study the case L = 2N, where 
the system is expected to be close to the transition be- 
tween a two- and a three-vortex state. Our preliminary 
findings indicate that the aforementioned mathematical 
structures may greatly simplify the CF study of this and 
other few-vortex states. 

The paper is organized as follows: We start by re- 
viewing the construction of CF trial wave functions in 
the context of rotating hard-core bosons, illustrated by 
some simple examples. We then present evidence that 
the overlap of the single vortex CF wave function with 
the exact one converges towards 1 in the large N limit. 
Since we believe that (modified versions of) our approach 
may provide a useful tool for studying other CF wave- 
functions with a high number of derivatives, our method 
of handling the derivatives is explained in some detail. In 
the following section, we discuss the region L < N. The 
CF wave functions are, in general, not exact here but do 
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reproduce the exact wave function for L = 0, 2, 3 which, 
again, suggests some interesting mathematical identities. 
FinaUy, we summarize our work. 



CF DESCRIPTION OF ROTATING BOSE 
SYSTEMS 

Consider a system of N spinless bosons with mass m 
in a harmonic trap of strength w, rotating with angular 
frequency and interacting via a short-range (delta func- 
tion) potential Hi. In a rotating frame the Hamiltonian 
can be written as 



N 



i=l 



2m 2 



ni^ + Hi 



(1) 



where Lz denotes the angular momentum around the ro- 
tation axis. Reshuffling a few terms, this is straightfor- 
wardly rewritten as 



N 



J-i^p.-AY + HUz. 



{uj - n)Lz + Hi{2) 



with A = muj{—y,x), Hho{z) denoting the z-part of the 
harmonic oscillator Hamiltonian and || denoting the pla- 
nar (x, y) part of the Hamiltonian. This is how the for- 
mal link to the quantum Hall system comes about: We 
see that the planar part of H takes the form of particles 
moving in an effective " magnetic" field B^f / = V x A = 
2mujz. Now, the interaction is assumed to be weak in the 
sense that it does not mix different harmonic oscillator 
levels. We will be interested, for a given total angular 
momentum, only in the lowest many-body states (the 
"yrast" band). In this limit, the model may be rewritten 
as a lowest Landau level (LLL) problem in the effective 
"magnetic" field B^ff = 2muj (and of course, = for 
the harmonic oscillator in the z-direction). The Hamil- 
tonian then takes the form 



H= {uj-VL)L + g^5^{vi-Yj) 



(3) 



(h = 1) where we now use L to denote the total angular 
momentum, L = li = L^. 

The single particle states spanning our Hilbert space 
(the lowest Landau level) are 



1 



110,1 = 



V2^+%n 



(4) 



where z = \/2muj{x + iy) are complex coordinates denot- 
ing the particle positions in the plane, and I is the an- 
gular momentum of the state. A general bosonic many- 
body wave function ip{zi, ...zn) will then be a homoge- 
neous, symmetric polynomial in the z^is, times the ex- 
ponential factor ea;p(— |zip/4) (which will be sup- 
pressed throughout this paper for simplicity). The degree 



of the polynomial gives the total angular momentum of 
the state. A special class of such wave functions are the 
so-called composite f ermion (CF) wave functions. They 
were first introduced by Jain 1C/| and have been extremely 
successful in describing FQH states, quantum dots in 
high magnetic fields and, more recently, h ighl y ro- 
tational states of Bose condensates |3 0, \m UJ. In 
quantum Hall physics, the main idea of this construction 
is, roughly speaking, to attach an even number of flux 
quanta to each electron, thus mapping them into weakly 
interacting composite fermions which can be thought of as 
moving in a reduced magnetic field. Technically, " attach- 
ing a flux quantum" means multiplying the wave function 
by a Jastrow factor. 



i<j 



(5) 



We see that the Jastrow factor has the effect of keep- 
ing the particles apart - it goes to zero if any two co- 
ordinates Zi and Zj approach each other. Therefore, it 
" takes care of much of the repulsive interaction between 
the particles. In the simplest approach, the so-called 
non-interacting composite fcrmion (NICF) approach, the 
wave function is then simply constructed as a Slater de- 
terminant of (non-interacting) composite fermions in the 
reduced magnetic field, times an even power of Jastrow 
factors. In the case of bosons, whose wave function has to 
be symmetric rather than antisymmetric, one instead ab- 
sorbs an odd power of flux quanta, mapping the bosons to 
weakly interacting composite fermions. In other words, 
bosonic trial wave functions with angular momentum L 
are constructed as non-interacting fermionic wave func- 
tions with angular momentum L — mN{N — l)/2, mul- 
tiplied by an odd number m of Jastrow factors (to = 1 
throughout this paper), and projected onto the LLL, 



(6) 



Here, fs denotes a Slater determinant consisting of 
single-particle wave functions rini{z,z) cx z^L\-^{zz/2) 
where n is the (CF) Landau level index (/ > —n) and 
LJj a generalized Laguerre polynomial. The LLL pro- 
jection V amounts to the replacement Zi 2d/dzi in 
the polynomial part of the wave function - the recipe is 
to replace all z:s with derivatives in the final polynomial, 
after multiplying out the Slater determinant and the Jas- 
trow factors and moving all z:s to the left. It has been 
shown. IQi] that with this projection method, the single- 
particle wave functions in the CF Slater determinant may 
be written as 



(7) 



with all derivatives acting only to the right. As this 
method tends to get computationally heavy in numerical 
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calculations with many particles and a large number of 
derivatives, slightly different methods of obtaining LLL 
wavefunctions have been employed in most of the CF 
literature[l3- These, too, are often referred to as projec- 
tion. Nevertheless, in this paper, "projection" will refer 
to the above brute force procedure. 

Before moving on to low angular momenta, let us illus- 
trate the method on two simple and well-known examples 
in the QH regime: First, consider the case L ~ N{N —1). 
Taking m — 1, the Slater determinant fs has to con- 
tribute an angular momentum N(N — l)/2 and is given 
by putting all CFs into the lowest CF Landau level, from 
l = Otol^ N -1, 
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EE - Z,). (8) 



i<j 



We immediately see from Eq. that the full wave func- 
tion is simply the bosonic Laughlin wave function 



tpL = Y[{zi ~ Zjf 



(9) 



with angular momentum L = N{N — 1). This is the ex- 
act ground state for the delta function interaction. Next, 
consider the angular momentum N{N — 1) — N, corre- 
sponding to a bosonic quantum Hall " quasiparticle" (as 
opposed to quasihole). In order to decrease the angular 
momentum by N as compared to the Laughlin state, we 
need to move one CF to the second Landau level, i.e. 
construct the Slater determinant 



fs = 



Zl 


Z2 


■ ZN 


1 


1 . 


1 


Zl 


Z2 


• Zpf 


Zl 


Z2 


Zn 


N-2 

Zl 


N-2 
Z2 


N-2 
■ Z]y 



(10) 



This gives the full trial wave function (again, apart from 
the exponential factor) 

N N 

'^gp = '^i^^y^t Yl i^k-zi) (z™ - z„) 

i=l k<hk,l=jii m<n 

^ 1 

EE;-3^n(-^--^)~'^^ (11) 

i—l j^i * k^i 



with ■0L denoting the Laughlin state Q • This wave func- 
tion has very high overlap with the exact one {e.g., 99.7 
% for 4 bosons [ll|). Its fermionic counterpart has been 
proven to capture correctly both the fractional charge 
and the anyonic statistics of the QH a uasielectron [lil lia | . 
and the same is expected to be the case for this bosonic 
version. 



Trial wave functions for other yrast states are con- 
structed in similar ways. The lower the angular momen- 
tum, the larger the number of derivatives. 



SINGLE VORTEX, L = N 

In the examples above, we considered angular mo- 
menta L ^ N"^ . Let us now turn to the regime L ^ N 
which, for large systems, corresponds to much smaller 
angular momenta. Note that since the Jastrow factor in 
Eq.(jHl has itself an angular momentum of raN{N — l)/2, 
we need to act with 0{N'^) derivatives in order to get 
down to L ^ N. One would expect that these deriva- 
tives acting on the Jastrow factor destroy most of the 
good correlations which arc at the very heart of the CF 
construction, which would thus fail in this regime. Sur- 
prisingly, however, we shall see that this is not the case. 
As shown in 111 , the Slater determinant for the ground 
state trial wave function at L = iV is constructed by oc- 
cupying the single-particle states ^n,— n 

for n = N ~ 2 

through n — 0, and the state ?/o,ii 
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(12) 



where, again, the z:s are to be replaced by derivatives 
acting on the Jastrow factor. Thus, the trial wave func- 
tion may be written as |lll | 



N 



V'(L=iV) = E(-l)"^" n idk~di)Y[iz.,~z,)il3) 

n=l k<l:k,l^n i<j 

On the other hand, it has been 

shown[3,ElE3 that the 
exact wave function for L = N , the single vortex, is given 
by 



N 



:N) 



(14) 



where Z = Zi/N is the center of mass. It was shown 
in Ref. [§| that, surprisingly, the overlap between ifT^ 
and 1)14(1 appeared to increase with increasing particle 
number (but only results up to = 10 were available at 
the time) . Here we have studied this point systematically 
and computed the overlaps for up to 43 particles. The 
mathematical approach which enabled us to handle such 
large numbers of derivatives, will be described shortly. 
First, however, let us state the result: The overlap be- 
tween the CF trial wave function J^l!^) and the exact one 
\lJi\j converges towards 1 in the large N limit, see fig. ^ 
The overlap equals 98.64% aheady at = 5 and 99.47% 
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for N = 10. For iV = 43 it equals 99.90%. We have 
checked that this is not a trivial result in the sense that 
the two wave functions might simply share the same lead- 
ing term and all other terms become irrelevant for large 
N. Rather, we found that the leading term of the two 
wave functions, Oi-^j' only corresponds to a weight of 
about 50% of the total wave function, and that the sub- 
leading terms of the two wave functions converge towards 
each other as well. This suggests that the functions ifT^ 
and (|14|l are identically equal to each other in the limit 
TV — > oo. This is a non-trivial mathematical identity 
which we have so far not been able to prove analytically; 
however, we noticed by simple curve-fitting, that the de- 
viation of the overlap from unity decreases as ~ 
We will discuss other, similar identities in the section on 
L <N. 



^.»******************* 
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FIG. 1: The overlap between the L = N CF trial wave func- 
tion 1131 1 and the exact one 1141 1 as function of the number of 
particles. 



Combinatorics and numerical method 



with 77 denoting all permutations over the particle in- 
dices (1,2,3,4), and e(7r) = ±1 for even/odd permuta- 
tions. Similarly, we rewrite the Jastrow factor as 



1111 

Zi Z2 Z3 Z4 

2 2 2 2 

Zi Z2 Z4 

Z3 ~3 ~3 -^3 
1 ^2 23 Z4 



— ^^'^i'^) ^ a (1)^1 (2) ^1(3)^1 {4,)- (17) 



Now, the first crucial observation is that an equivalent 
way of writing the permutation sum (|17|l is by fixing the 
particle indices and summing over all permutations of the 
exponents instead. 



J = Ve(CT) z° 



-^Ca)-! <t(3)-1 <t(4)-1 



(18) 



Moreover, since reordering the terms will leave the sum 
unchanged, we may apply any permutation tt to all the 
indices as long as we make sure to change the signs ac- 
cordingly: 



cr(l)-l ct(2)-1 cr(3)~l (t(4)-1 



7r(l) 7r(2) 



^t(3) ^7.(4) 



(19) 



Noting that e(7r)e(7r(T) = e{a), we can now combine the 
Slater determinant fs and the Jastrow factor J to give 
the following expression for the (unnormalized) single 
vortex wave function, 



(4) 



'^(2)-l -^CS)-! a(4)-l 



t(1) -7,(2) -7r(3) 



(20) 



We note that the outer sum is simply a (symmetric) sum 
over all permutations tt of the particles indices. There- 
fore, it is sufficient to study one term, denoted as fu, 
in this permutation sum and perform the symmetriza- 
tion only at the very end. Choosing this term to be 
^= (1,2,3,4), we get 



fu 



E 



eia)z, d, dl X zl'^'^-'z'^.^'^-'zf^-'zl^'^-' 



For the sake of clarity, we will illustrate our approach 
to the single vortex wave function for the simple case N = 
4. The generalization to arbitrary N is straightforward, 
and the general result will be stated at the end. After 
projection, the CF Slater determinant H12|l takes the form 



fs 



Zl Z2 Z-i Z4 
1111 

9i di 94 
dl dl dl dl 



(15) 



This determinant may be rewritten as a permutation 
sum. 



/s = ^e(^)zi(i)Z°(2)5i(3)9^ 



(4) 



(16) 



^e(cr) zl 



(a(3) - l)z: 



(t(3)-2 



X (a(4)-l)(a(4)-2)zl(^)-^ 



(21) 



The next important observation is that the above sum 
may again be expressed as a determinant. 



fu^ 



Zl 1 
zl Zl 




1 

2z3 
2 
3 





2 

Z2 3z| 6z4 



(22) 



Thus, to summarize, we have expressed both the Slater 
determinant and the Jastrow factor in lO as permutation 
sums, combined these, performed all derivatives, and cast 
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the final expression in the form of a new determinant. 
The full wave function '4'iL=N) is obtained by adding all 
permutation of the particle indices in (|22|l . We have ob- 
tained an explicit form of the CF polynomial, without 
derivatives, which may then be used in numerical calcu- 
lations, to compute overlaps etc. 

The generalization of H22|l to arbitrary N is straight- 
forward; one finds that 



of basis we can therefore convert 4'('l=n) *° orthogo- 
nal symmetrized basis, enabling us to compute the over- 
laps. 

The method described in this section may be straight- 
forwardly modified to other values of L. In the next 
section we illustrate how it may be used to shed light on 
some yrast states at even lower angular momenta. 



Zl 

4 



1 

Z2 



1 



(23) 



A particularly appealing property of this matrix is its 
near-triangular form, which greatly reduces the calcula- 
tion cost. 

Before ending this section, let us briefly outline how 
we used this result to compute overlaps between the CF 
wave function and the exact one l|14|l. The polynomial 
part of '0(L=Af) may be regarded as a vector in the space of 
symmetric polynomials. A convenient basis for this space 
consists of the symmetrized sums of individual polyno- 
mial terms. For example, the symmetrization of ziZ2z| 
is 



S[ziZ2z'l] = 2ziZ2z| + 2Z1Z2Z4 + 2ziZ3z| 



(24) 



and may be regarded as one basis vector. (The factor 
2 arises from the fact that both eg. ziZ2z'^ and 222123 
occur in the sum.) Since ■(/'(L=Af) = S[fu] {S denoting 
symmetrization), it is therefore straightforward to deter- 
mine the representation of ip^L^N) in tii6 symmetric basis 
by symmetrizing each term of fjj individually. 

Since overlap integrals between polynomial terms with 
non-matching powers are zero, our symmetric basis is or- 
thogonal. This means that if both ip(^L=N) and V'(i=jv) 
are normalized and represented in a normalized symmet- 
ric basis, then the overlap may be computed as a simple 
dot product. This also implies that all calculated overlap 
values are exact (up to numerical precision). To express 
''P{L=N) ^^^^ basis, let us rephrase the product H14I) as 
the sum 



^ 1 . 



(25) 



where Si = 3^(21, . . . ,zm) are the fundamental symmetric 
polynomials of degree i , 



i{zi, . . . , Zn) = S [ziZ2...Zi 



(26) 



The products sat, siS(Ar_i), ■ ■ ■ ,3^ appearing in this sum 
form another, non-orthogonal, basis. Through a change 



BELOW THE SINGLE VORTEX, L <N 

Exact ground state wave functions arc known not only 
for the single vortex (|14() but for all angular momentum 
states 2 < L < N. As was shown some years ago [l^ll5j . 
they are simply given by fundamental symmetric poly- 
nomials SLizi) where Zi = Zi — Z and Z is again the 
center-of-mass coordinate: 



E 

Pi<P2<-<Pl 



{zpi - Z)izp^ - Z)--- [zp^ - Z).{21) 



For example, 'ipL=2 — S [{zi — Z){z2 — Z)], with S denot- 
ing symmetrization over all particle coordinates. Note 
that these states are translation invariant (TI), i.e. in- 
variant under a simultaneous constant shift Zi Zi + a 
of all the coordinates^,!!^!. We see that the single vor- 
tex (|14|l is merely a special case of this series of wave 
functions. (At L = the polynomial part of the exact 
ground state is trivially a constant as all particles are in 
the I — state, while the only possible way of construct- 
ing the ground state at L = 1 is ^^^L=l — J2i -^ii which 
corresponds to a center of mass excitation of the L ~ 
state.) 

As discussed in Refs. 0,113, the space of TI polyno- 
mials for 2 < i < is spanned by the basis states 



(28) 



where L — X]n=2'^^n- th.us see that the states (|27f) 
are special cases of such basis states, with — 1, all 
other fc„ = 0. In anticipation of the discussion below, 
note that for L = 2,3, they are the only possible basis 
states, and thus it is obvious that (|27() has to be exact 
for L = 2, 3. (For larger L this is no longer the case; for 
example, the basis for L — 4 consists of |20000...) as well 
as 100100...).) 

Given these results, let us discuss some mathematical 
peculiarities of the CF construction in the regime L < N. 
First of all, the CF construction immediately reproduces 
the exact wave function for L = 0: In order to cancel 
all the angular momentum of the Jastrow factor, the CF 
Slater determinant with lowest possible (quasi) Landau 
level energy is uniquely given by occupying the CF Lan- 
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dau level states ?7„,-„ for n = through N — 1, 



fs = 
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1 . 
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Zl 














Z2 


Zn 




-N-1 


-N-1 




z-i 


■ Zj\[ 



= Y[iz.- (29) 



After projection, this will give a Jastrow factor of deriva- 
tives acting on the corresponding Jastrow factor of z:s, 
thus resulting in a constant. This can be made explicit 
using the notation of the previous section. Proceeding 
exactly as in equations (|16l) - (|22|l . again taking = 4, 
we find that Vl=o = fui^l), 7^(2), 7r(3), 7r(4)) with 



fu = J2e{a)d',dldidl X zl^'^z^,^'^zf^zl^'-^ (30) 
= e{a)zl^'^ ■ o{2)zl^^-^ ■ a(3)(a(3) - V)zl^'^-^ 



a(4)(a(4)-l)(a(4)-2)z. 



<7(4)-3 



(31) 



where a now denotes all permutations over (0,1,2,3). 
This is identified with the determinant 



1 





1 
zl 2z2 2 
zl 3z| 623 6 



12 



(32) 



which obviously gives the correct wave function (up to 
normalization) for L = (all bosons in the I = state). 

The CF construction reproduces the exact ground state 
wave functions for L = 2 and 3 as well. In contrast to 
the case L — N, these identities are not asymptotic - 
they are exact for all N. This is not surprising, for the 
following reason [l]|: For these angular momenta, one 
can construct compact CF wave functions which are thus 
known to be translation invariant (TI) 21]. On the other 
hand, as discussed above, for these lowest angular mo- 
menta there only exists one basis state in the space of 
TI polynomials. Therefore, the CF compact state has to 
correspond to this basis state, which is the exact wave 
function. 

What is somewhat surprising, however, is how this 
comes about mathematically. Consider again the L = 
state discussed above. From this we can construct com- 
pact L = 2 states by moving any of the particles (except 
the one at (0, 0)) from the CF Landau level state (n, —n) 
to (n — 1, —n + 2). There are thus — 1 different ways 
of constructing the L = 2 CF Slater determinant which 
are all degenerate in CF kinetic energy. In mathemati- 
cal terms, this corresponds to replacing one of the states 
Vn,-n = 9" by (see (U))). It thus looks as if there 

are iV — 1 different CF candidate wave functions for the 
L = 2 state, while according to the above arguments, we 
would expect the CF wave function to be unique and ex- 
act! The solution to this apparent paradox is that while 



the N — 1 possible Slater determinants are certainly not 
equal to each other, they all result in the same polynomial 
after acting on the Jastrow factor. To see how this comes 
about, consider again the case of four particles. The first 
possible Slater determinant is 



Zl Z2 Z3 Z4 
1111 

dl dl dl dl 
dl 82 (91 (91 



(33) 



resulting in 

fu = Y.<a)z,zldldl X zl^'^z^^'^f^zl^'^ (34) 
= j2^{o)zl^"'^'.zl^'^.ama{i)-l)zf^-' 



■ a(4)(a(4)-l)(a(4)-2)z4 
which equals the determinant 



o-(4)-3 



fu 



Zl 1 
zl Z2 







2 

6Z3 6 



12(ziZ2 - zl). 



(35) 



(36) 



Proceeding in the same way for the other two possibilities 
gives fu = 12(2! - Z2Z3) and fu = 12(z| - 2:324), respec- 
tively. Thus, since the full wave function is constructed 
from fu by summing over all permutations of the indices 
(1,2,3,4), the final, symmetrized sum will be the same for 
all these constructions. 

Exactly the same thing happens for L = 3. There 
are N — 2 ways of constructing the Slater determinant, 
obtained from the L = one by letting any 9" — ^ 2(}"~^. 
Again, these different constructions simply correspond 
to picking out different parts of the final, symmetrized 
polynomial. 

Of course, similar things may happen at other angu- 
lar momenta as well. For example, there are two ways 
of constructing the Slater determinant for L = — 1. 
Again, these result in the same polynomial, though in 
this case the CF construction does not give the exact 
wave function even in the large N limit. 

The fact that the wave functions resulting from the CF 
construction for L = 2, 3 are equal to the very simple and 
compact expressions in (|27f) . again provides non-trivial 
and potentially useful mathematical identities. 



SUMMARY AND OUTLOOK 

To summarize, we have shown that the CF construc- 
tion for bosons at low angular momenta may be taken 
to far higher particle numbers than previously by han- 
dling the derivatives of the CF Slater determinant in a 
new way. In particular, we used this to show that for 
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the single vortex, L = N, the CF construction becomes 
exact in the high N hmit, leading to a non-trivial and po- 
tentially useful mathematical identity. For other L, the 
CF construction does, in general, not produce the exact 
wave functions. We have argued why it is expected to 
give the exact ground state for L = 2,3 (all N). How- 
ever, work remains to get a better understanding of what 
makes L = N special. Moreover, we used our methods 
to show examples of other mathematical peculiarities of 
the CF construction in the regime L < N . Our hope is 
that the results of this paper may provide new insight 
into properties of CF wave functions in the few-vortex 
regime beyond L = N, where no exact, analytical wave 
functions are known. For example, a preliminary study 
of the case L — 2N suggests that while the number of 
possible CF Slater determinants increases linearly with 
N, these may reduce to just two different wave functions; 
finding the ground state as a linear combination of these 
would then just be a one-parameter problem. There is 
thus hope that one may construct rather simple analyti- 
cal trial wave functions in the regime L ^ N , which may 
directly reveal the vortex structures of these states. The 
present paper indicates that in this regime, the CF con- 
struction may do far better than previously assumed in 
the literature. 
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